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We study the hadron-quark phase transition in the interior of neutron stars (NS). For the hadronic 
sector, we use a microscopic equation of state (EOS) involving nucleons and hyperons derived within 
the Brueckner-Bethe-Goldstone many-body theory, with realistic two-body and three-body forces. 
For the description of quark matter, we employ both the MIT bag model with a density dependent 
bag constant, and the color dielectric model. We calculate the structure of NS interiors with the 
EOS comprising both phases, and we find that the NS maximum masses are never larger than 1.7 
solar masses, no matter the model chosen for describing the pure quark phase. 
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PACS numbers: 26.60.+C, 21.65.+f, 24.10.Cn, 97.60.Jd 



I. INTRODUCTION 

The appearence of quark matter in the interior of mas- 
sive neutron stars (NS) is one of the main issues in the 
physics of these compact objects. Calculations of NS 
structure, based on a microscopic nucleonic equation of 
state (EOS), indicate that for the heaviest NS, close to 
the maximum mass (about two solar masses), the central 
particle density reaches values larger than 1/fm 3 . In this 
density range the nucleon cores (dimension « 0.5 fm) 
start to touch each other, and it is hard to imagine that 
only nucleonic degrees of freedom can play a role. On 
the contrary, it can be expected that even before reach- 
ing these density values, the nucleons start to loose their 
identity, and quark degrees of freedom are excited at a 
macroscopic level. 

Unfortunately it is not straightforward to predict the 
relevance of quark degrees of freedom in the interior of 
NS for the different physical observables, like cooling evo- 
lution, glitch characteristics, neutrino emissivity, and so 
on. In fact, the other NS components can mask the ef- 
fects coming directly from quark matter. In some cases 
the properties of quark and nucleonic matter are not very 
different, and a clear observational signal of the presence 
of the deconfined phase inside a NS is indeed hard to 
find. 

The value of the maximum mass of NS is probably 
one of the physical quantities that are most sensitive to 
the presence of quark matter in NS. If the quark mat- 
ter EOS is quite soft, the quark component is expected 
to appear in NS and to affect appreciably the maximum 
mass value. In this case the maximum mass is expected 
to be slightly larger than the observational limit (1.44 
solar masses of the so-called Taylor pulsar pj). The ob- 
servation of a large NS mass (larger than 2 solar masses) 
would imply, on the contrary, that the EOS of NS mat- 
ter is stiff enough to keep the maximum mass at this 
large values. Purely nucleonic EOS are able to accomo- 
date masses comparable with these large values 
Since the presence of non-nucleonic degrees of freedom, 
like hyperons and quarks, tends usually to soften consid- 
erably the EOS with respect to purely nucleonic matter, 
thus lowering the mass value, their appearence would be 



incompatible with observations. The large value of the 
mass could then be explained only if both hyperonic and 
quark matter EOS are stiffer than expected. In particu- 
lar, the quark EOS should be assumed to be stiff enough 
to render the deconfined phase energetically disfavoured. 

In this paper we will discuss this issue in detail. Un- 
fortunately, while the microscopic theory of the nucle- 
onic EOS has reached a high degree of sophistication, 
the quark matter EOS is poorly known at zero temper- 
ature and at the high baryonic density appropriate for 
NS. One has, therefore, to rely on models of quark mat- 
ter, which contain a high degree of uncertainity. The 
best one can do is to compare the predictions of differ- 
ent models and to estimate the uncertainty of the results 
for the NS matter as well as for the NS structure and 
mass. In this paper we will use a definite nucleonic EOS, 
which has been developed on the basis of nuclear mat- 
ter many-body theory, and two different models for the 
quark EOS, and compare the results. Confrontation with 
previous calculations shall also be discussed. 

The paper is organized as follows. In section [H] we 
review the determination of the baryonic EOS compris- 
ing nucleons and hyperons in the Brucckncr-Hartree-Fock 
approach. Section IIIII concerns the quark matter EOS 
according to the MIT bag model and the color dielec- 
tric model (CDM). In section llVl we present the results 
regarding neutron star structure combining the baryonic 
and quark matter EOS for beta-stable nuclear matter. 
Section IV! contains our conclusions. 



II. BRUECKNER THEORY 
A. EOS of nuclear matter 

The Brueckner-Bethe-Goldstone (BBG) theory is 
based on a linked cluster expansion of the energy per 
nucleon of nuclear matter (see Ref . [5j , chapter 1 and ref- 
erences therein). The basic ingredient in this many-body 
approach is the Brueckner reaction matrix G, which is 
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the solution of the Bethe-Goldstone equation 
^ \k a h) Q (k a kb\ 



G[p; uj] = v - 



G[p;c 



(1) 



where v is the bare nucleon-nucleon (NN) interaction, p 
is the nucleon number density, and to the starting energy. 
The single-particle energy e(k) (assuming Ti—1 here and 
throughout the paper), 



e(fc) = e(fc;p) = — + U(k;p), 
2m 



(2) 



and the Pauli operator Q determine the propagation 
of intermediate baryon pairs. The Brueckner-Hartree- 
Fock (BHF) approximation for the single-particle poten- 
tial U(k; p) using the continuous choice is 

U(k;p)=Re (kk'\G[p;e(k) + e(k')]\kk') a , (3) 

k'<k F 

where the subscript "a" indicates antisymmctrization of 
the matrix element. Due to the occurrence of U{k) in 
Eq. J5J, they constitute a coupled system that has to 
be solved in a self-consistent manner for several Fermi 
momenta of the particles involved. In the BHF approxi- 
mation the energy per nucleon is 



E _ 3 fc| 1 
A~52^ + 2p~ ^ 



(kk'\G[p;e(k)+e(k')]\kk') a . (4) 



In this scheme, the only input quantity we need is the 
bare NN interaction v in the Bethe-Goldstone equation 
In this sense the BBG approach can be considered 
as a microscopic one. The nuclear EOS can be calcu- 
lated with good accuracy in the Brueckner two hole-line 
approximation with the continuous choice for the single- 
particle potential, since the results in this scheme are 
quite close to the calculations which include also the three 
hole-line contribution [6|. In the calculations reported 
here, we have used the Argonne vi% potential J7j as the 
two-nucleon interaction. 

However, it is commonly known that nonrelativistic 
calculations, based on purely two-body interactions, fail 
to reproduce the correct saturation point of symmetric 
nuclear matter, and three-body forces (TBF) among nu- 
cleons are needed to correct this deficiency. In this work 
the so-called Urbana model will be used, which consists 
of an attractive term due to two-pion exchange with ex- 
citation of an intermediate A resonance, and a repulsive 
phenomenological central term [8j. We introduced the 
same Urbana three-nucleon model within the BHF ap- 
proach (for more details see Ref. 0). In our approach the 
TBF is reduced to a density dependent two-body force by 
averaging over the position of the third particle, assum- 
ing that the probability of having two particles at a given 
distance is reduced according to the two-body correlation 
function determined self-consistently. The corresponding 
nuclear matter EOS fulfills several requirements, namely 



(i) it reproduces the correct nuclear matter saturation 
point, (ii) the incompressibility is compatible with the 
values extracted from phenomenology, (iii) the symme- 
try energy is compatible with nuclear phenomenology, 
(iv) the causality condition is always fulfilled. 

In order to study the structure of neutron stars, we 
have to calculate the composition and the EOS of cold, 
neutrino-free, catalyzed matter. We require that the neu- 
tron star contains charge neutral matter consisting of 
neutrons, protons, and leptons (e~, p~) in beta equi- 
librium, and compute the EOS for charge neutral and 
beta-stable matter in the following standard way 
The Brueckner calculation yields the energy density of 
lepton/baryon matter as a function of the different par- 
tial densities, 
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(5) 



where we have used ultrarelativistic and nonrelativistic 
approximations for the energy densities of electrons and 
muons, respectively In practice, it is sufficient to com- 
pute only the binding energy of symmetric nuclear matter 
and pure neutron matter, since within the BHF approach 
it has been verified 10, 11] that a parabolic approxima- 
tion for the binding energy of nuclear matter with ar- 
bitrary proton fraction x = p p / p, p = p n + p p , is well 
fulfilled, 



0.5) + (1 - 2.x) 2 £; sym (p) 



(6) 



A A' 

where the symmetry energy E sym can be expressed in 
terms of the difference of the energy per particle between 
pure neutron (x = 0) and symmetric (x — 0.5) matter: 



E sym (p) - - 



1 d(E/A) 



( P ,0)«^(p,0)-|(p,0.5). (7) 



4 dx vr ' ' A 
Knowing the energy density Eq. (|SJ| , the various chem- 
ical potentials (of the species i — n,p,e,p) can be com- 
puted straightforwardly, 

and the equations for beta-equilibrium, 

Pi = bip n - q. t p e , (9) 

(hi and qi denoting baryon number and charge of species 
i) and charge neutrality, 



E 



PiQ, 



= 0. 



(10) 



allow one to determine the equilibrium composition pi(p) 
at given baryon density p and finally the EOS, 

d e(pi(p)) de 



P{P) = P' 



dp p 



P-r -e = ppn-e 
dp 



(11) 
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B. Hyperons in nuclear matter 



uous choice are given by 



While at moderate densities p ss p the matter in- 
side a neutron star consists only of nucleons and lep- 
tons, at higher densities several other species of particles 
may appear due to the fast rise of the baryon chemi- 
cal potentials with density. Among these new particles 
are strange baryons, namely, the A, E, and 5 hyperons. 
Due to its negative charge, the S~ hyperon is the first 
strange baryon expected to appear with increasing den- 
sity in the reaction n + n — ► p + E~ , in spite of its sub- 
stantially larger mass compared to the neutral A hyperon 
(Af E - = 1197 MeV,M A = 1116 MeV). Other species 
might appear in stellar matter, like A isobars along with 
pion and kaon condensates. It is therefore mandatory to 
generalize the study of the nuclear EOS with the inclu- 
sion of the possible hadrons, other than nucleons, which 
can spontaneously appear in the inner part of a neutron 
star, just because their appearance is able to lower the 
ground state energy of the dense nuclear matter phase. 

As we have pointed out in the previous section, the nu- 
clear EOS can be calculated with good accuracy in the 
Brueckner two hole-line approximation with the contin- 
uous choice for the single-particle potential, since the re- 
sults in this scheme are quite close to the full convergent 
calculations which include also the three hole-line contri- 
bution. It is then natural to include the hyperon degrees 
of freedom within the same approximation to calculate 
the nuclear EOS needed to describe the neutron star in- 
terior. To this purpose, one requires in principle nucleon- 
hyperon (NY) and hyperon-hyperon (YY) potentials. In 
our work we use the Nijmegen soft-core NY potential 
[T3 | that is well adapted to the available experimental 
NY scattering data. Unfortunately, up to date no YY 
scattering data and therefore no reliable YY potentials 
are available. We therefore neglect these interactions in 
our calculations, which is supposedly justified, as long as 
the hyperonic partial densities remain limited. Also, for 
the following calculations the vis NN potential together 
with the phenomenological TBF introduced previously, 
are used. 

With the NN and NY potentials, the various G matri- 
ces are evaluated by solving numerically the Brueckner 
equation, which can be written in operatorial form as 

B0 



G ab [W] = V ab + J2 E V -c\PP') 
c p,p f 



Qc 



W - E c 



-( PP '\G cb [W}, 



(12) 

where the indices a, 6, c indicate pairs of baryons and the 
Pauli operator Q and energy E determine the propaga- 
tion of intermediate baryon pairs. In a given nucleon- 
hyperon channel c = (NY) one has, for example, 



E(NY) = m N + m Y + X 



k'y 



2uin 2m Y 



+ U N (k N ) + U Y (k Y ). 
The hyperon single-particle potentials within the contin- 



Uy(k) = Re ]T E 

N = n 'Pk'<ki N} 



x(kk'\G {NYKNY) [E {NY) (k,k')] kk') (14) 



and similar expressions of the form 

r\(/,)= >. (V '"•■) + E ug\k) (15) 



E <>(*) 

N'=n,p 



K=£-,A 



apply to the nucleon single-particle potentials. The nu- 
cleons feel therefore direct effects of the other nucleons 
as well as of the hyperons in the environment, whereas 
for the hyperons there are only nucleonic contributions, 
because of the missing hyperon-hyperon potentials. The 
equations (|12f) -H15 [) define the BHF scheme with the con- 
tinuous choice of the single-particle energies. In contrast 
to the standard purely nucleonic calculation there is now 
an additional coupled channel structure, which renders a 
self-consistent calculation quite time-consuming. 

Once the different single-particle potentials are known, 
the total nonrelativistic baryonic energy density, e, can 
be evaluated: 



E 



k * dk k 2 



i— n,p,5H — ,A 

ejv« + E 



k 2 1 „ 
2^+2^ 



(16) 
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K=S-,A 
k 2 



! * dk k 2 



my + 



2m Y 



+ U Y n \k)+Uy V) (k) 



(17) 



where eatat is the nucleonic part of the energy density, 
Eq. Using for example an effective mass approxima- 
tion for the hyperon single-particle potentials, one could 
write the last term due to the nucleon-hyperon interac- 
tion as 



E 

y=E-, 



Py 



m Y + U Y (0) 



1 (3^py) 5 / 3 



2ml 



5tt 2 



(18) 

which should be added to Eq. JSJ. 

The knowledge of the energy density allows one then 
to compute EOS and neutron star structure as de- 
scribed before, now making allowance for the species 
i = n, p, Yi~ , A, e~ , pT . The main physical features of the 
nuclear EOS which determine the resulting compositions 
are essentially the symmetry energy of the nucleon part 
of the EOS and the hyperon single-particle potentials in- 
side nuclear matter. Since at low enough density the nu- 
cleon matter is quite asymmetric, the small percentage 
of protons feel a deep single-particle potential, and there- 
fore it is energetically convenient to create aE" hyperon, 
since then a neutron can be converted into a proton. The 
depth of the proton potential is mainly determined by the 
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nuclear matter symmetry energy. Furthermore, the po- 
tentials felt by the hyperons can shift substantially the 
threshold density at which each hyperon sets in. 

We have found rather low hyperon onset densities of 
about 2 to 3 times normal nuclear matter density for 
the appearance of the S~ and A hyperons ^j. (Other 
hyperons do not appear in the matter). Moreover, an 
almost equal percentage of nucleons and hyperons are 
present in the stellar core at high densities. The inclusion 
of hyperons produces an EOS which turns out to be much 
softer than the purely nucleonic case. The consequences 
for the structure of the neutron stars are dramatic. In 
fact the presence of hyperons leads to a maximum mass 
for neutron stars of less than 1.3 solar masses [Tojj . which 
is below the observational limit. 

This surprising result is due to the strong softening of 
the baryonic EOS when including hyperons as additional 
degrees of freedom. We do not expect substantial changes 
when introducing refinements of the theoretical frame- 
work, such as hyperon-hyperon potentials 14], relativis- 
tic corrections, etc. Three-body forces involving hyper- 
ons could produce a substantial stiffening of the baryonic 
EOS. Unfortunately they are essentially unkown, but can 
be expected to be weaker than in the non-strange sector. 
Another possibility that is able to produce larger maxi- 
mum masses, is the appearence of a transition to another 
phase of dense (quark) matter inside the star. This will 
be discussed in the following. 



III. QUARK PHASE 

The results obtained with a purely baryonic EOS call 
for an estimate of the effects due to the hypothetical 
presence of quark matter in the interior of the neutron 
star. Unfortunately, the current theoretical description 
of quark matter is burdened with large uncertainties, se- 
riously limiting the predictive power of any theoretical 
approach at high baryonic density. For the time being 
we can therefore only resort to phenomenological models 
for the quark matter EOS and try to constrain them as 
well as possible by the few experimental information on 
high density baryonic matter. 

One of these constraints is the phenomenological ob- 
servation that in heavy ion collisions at intermediate en- 
ergies (10 MeV/A < E/A < 200 MeV/A) no evidence 
for a transition to a quark-gluon plasma has been found. 
Indeed, all microscopic simulations, like BUU or QMD, 
that are able to reproduce a great variety of experimental 
data, do not need the introduction of such a transition 
|l5|. I n these simulations the calculated nucleon density 
can reach values which are at least 2 to 3 times larger 
than the saturation density po- One can, therefore, con- 
clude that symmetric or nearly symmetric nuclear matter 
at a few MeV of temperature does not exhibit any phase 
transition to deconfined matter up to this baryon density. 
It has to be noticed that the phase transition in symmet- 
ric matter can occur at a substantially different baryon 



density than in neutron star matter, where nuclear mat- 
ter is closer to neutron matter than to symmetric nuclear 
matter. 

This constraint coming from heavy-ion physics appears 
as an independent one, that should be fulfilled by any 
theory or model of deconfmement. Indeed, quark matter 
models can have, in some cases, serious difficulties to 
fulfill the constraint (the transition occuring at too low 
density), even if they produce "reasonable" results for 
neutron stars, where the transition does occur. We will 
in the following take this constraint in due consideration, 
and use an extended MIT bag model ^(| (including the 
possibility of a density dependent bag "constant") and 
the color dielectric model [l7| , both compatible with this 
condition. 

Furthermore, some theoretical interpretation of the 
heavy ion experiments performed at the CERN SPS 0] 
points to a possible phase transition at a critical den- 
sity p c ss 6po ss 1/fm 3 . This condition has to be taken 
with caution, since in ultrarelativistic heavy ion colli- 
sions, quark matter is probed essentially at zero baryonic 
density and high temperature, at variance with the case 
of NS matter, where high baryonic density and essentially 
zero temperature are present. 



A. The MIT bag model 

We first review briefly the description of the bulk prop- 
erties of uniform quark matter, deconfined from the (3- 
stable hadronic matter mentioned in the previous section, 
by using the MIT bag model ^(|. The thermodynamic 
potential of / = u,d,s quarks can be expressed as a 
sum of the kinetic term and the one-gluon-exchange term 
flfj I2H ] proportional to the QCD fine structure constant 



3m,j 
s 2tt 3 



Vf x f 



(2a3-3)+ln(s/ + V/) 
1 2 2 

yf x f - M x f + Vf) 



My/) 

(19) 



where m f and p, f are the / current quark mass and chem- 
ical potential, respectively, yf — fif jra f, xj — \Jy 2 f — 1, 

and CT ron = 313 MeV is the renormalization point. In this 
work we will consider massless u and d quarks (together 
with m s = 150 MeV), in which case the above expression 
reduces to 



n q = 



-SU-^)- ( q = u,d). (20) 



4tt 2 



The number density pj of / quarks is related to f2/ via 



PS 



dp f ' 



(21) 
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and the total energy density for the quark system is writ- 
ten as 



£mit(p ill Pdi Ps 



f 



(Q f + njpf) + B 



(22) 



where B is the energy density difference between the 
perturbative vacuum and the true vacuum, i.e., the 
bag "constant." In the original MIT bag model B sa 
55 MeVfar 3 is used, while values B w 210 MeVfm -3 
are estimated from lattice calculations [2l] . In this sense 
B can be considered as a free parameter. 

The composition of /3-stable quark matter is deter- 
mined by imposing the condition of equilibrium under 
weak interactions for the following processes 



u + e 
u + e" 
d 
s 

s + u 



d + L> e , 
S + V e , 

u + e~ + V e , 
u + e~ + V e , 
d + u . 



(23a) 
(23b) 
(23c) 
(23d) 
(23e) 



In neutrino- free matter (p,^ = pv c = 0), the above equa- 
tions imply for the chemical potentials 



Pd = Ps = Pu + Pe 



(24) 



As in baryonic matter, the relations for chemical equilib- 
rium must be supplemented with the charge neutrality 
condition and the total baryon number conservation, 



= g(2p„ -pd~ Ps) - Pe , 
P = g(/>« + Pd + Ps) , 



(25) 
(26) 



in order to determine the composition Pf(p) and the pres- 
sure of the quark phase, 



(27) 



It has been found |2J, |23j that within the MIT bag 
model (without color superconductivity) with a density 
independent bag constant B, the maximum mass of a 
NS cannot exceed a value of about 1.6 solar masses. In- 
deed, the maximum mass increases as the value of B 
decreases, but too small values of B are incompatible 
with a transition density p > (2, . . . , 3)po in symmetric 
nuclear matter, as demanded by heavy-ion collision phe- 
nomenology according to the preceeding discussion. In 
order to avoid this serious drawback of the model, one 
can introduce a density-dependent bag "constant" B(p), 
and this approach was followed in Refs. [2^] . This allows 
one to lower the value of B at large density, providing 
a stiffer quark matter EOS and increasing the value of 
the maximum mass, while at the same time still fulfill- 
ing the condition of no phase transition below p w 3po. 
The comparison of the predictions based on the MIT bag 



model and the CDM can be considered meaningful only 
if this constraint is maintained also in the CDM. 

In the following we present results based on the MIT 
model using a constant value of the bag constant, B = 
90 MeV/fm 3 , and a gaussian parametrization for the den- 
sity dependence, 



B(p)=B 00 + (B Q -B 00 )cxp 



X-) 



(28) 



with Boo = 50 MeVfrn" 3 , B = 400 MeVfm~ 3 , and 
(3 = 0.17. 



B. The color dielectric model 

Let us now consider the color dielectric model, which 
was originally introduced |17| as a confinement model of 
the nucleon (see the general reviews 0, |2jJ |2(|). In 
the CDM the nucleon is described as a soliton in which 
quarks are dynamically confined via the interaction with 
a scalar-isoscalar chiral singlet field, indicated in the fol- 
lowing as x, whose quanta correspond to glueballs or hy- 
brid mesons. Several closely related versions of the CDM 
exist in the literature. They have been widely employed 
in the baryon sector, to calculate the static p7 ll28ll29ll30| 
and dynamic properties 0, |U H3, 113 l35Tp6l o f the nu- 
cleon and to describe strange barvons |28l l37l l38l \3^ . 
The CDM has also been applied in the quark sector . to 
calculate the EOS of quark matter HO, EE |H lil lil 
and to study the stability of strange quark matter |45| . 
Applications of the CDM EOS for quark matter to the 
study of compact stars have been considered by Ghosh et 
al. |46J and by Drago and collaborators, who studied the 
structure of hybrid stars 0; HH an d the problem of su- 
pernova explosions |49| . and, more recently by Malheiro 
et al. |5Ct l5l| to study pure quark stars. 

In this work we use the chiral version of the CDM 
0, extended to include strange quarks 0,E3> which, as 
shown in Ref . |44| , describes reasonably well the nucleon, 
while giving, with the same set of parameter values, a 
meaningful equation of state for symmetric quark matter. 
The Lagrangian of the model reads 

c = E WffWi 

f—u,d,s 

+ Y] -T-^t (f + iT • 7T7 5 ) ipf - — ll> s ll)s 



f=u,d 



+ \{d»<jf+ l -(d^) 2 -U{<j^) 
+ \{d,x) 2 -V( X ). 
Here the potential V(x) nas the quadratic form 

v(x) = \m\ 2 , 



(29) 



(30) 
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while 



-ft) 



(31) 



is the usual "Mexican hat" potential. The character- 
istic feature of the CDM is the coupling of the quarks 
to an inverse power of the field Xj through which the 
quarks acquire density dependent effective masses m Ut d = 
— g u ,dG / f-nX an( i m s — 9s/X- Since x vanishes in the vac- 
uum, the quark masses diverge as the density p goes to 
zero, thus providing confinement. 

In Eq. lt2~9"|) the couplings are given by g Ut d = #(,^±£3) 
and g s = g{2f K - f n ), where f n = 93 MeV and f K = 113 
MeV are the pion and the kaon decay constants, respec- 
tively, and where £3 = f K ± — /r-o = —0.75 MeV. In 
Eq. gTJ we take m a = 1.2 GeV. 

The parameters of the model are thus g, determining 
the couplings <?/, and the mass M of the x field. At 
the mean field level, and with the form Eq. I|30(l of the 
potential, the only free parameter is actually the prod- 
uct G — \fg~M. In the following we will use the val- 
ues M — 1.7 GeV and g — 23 MeV (corres pon ding to 
G = 198 MeV), which were obtained in Ref. H3, by re- 
quiring that the model provides reasonable values for the 
average delta-nucleon mass and for the nucleon isoscalar 
radius. 

To describe the possible quark phase in the interior of a 
neutron star we consider a uniform system of plane wave 
u, d, and s quarks interacting with the fields x an( i <J - 
In the mean field approximation the latter are assumed 
to be constant, while the pion field vanishes. The energy 
density of the system is then given by 



ccdm = 7 Y Y V fe2 + m /( cr 'X) 2 

f=u,d,s k<k f F 

+ V( X ) + U(a,w = 0), (32) 

where 7 = 6 is the spin and color degeneracy factor and 
kp are the Fermi momenta of the quarks of flavor /. At 
fixed baryon density, two coupled equations for the fields 
X and a are obtained by minimizing the energy density 
£cdm: 



Ps( k F' m f) 



d-X f-^. 

f=u,d 

dU(*,0) x . , , 



9f°_ 

Ux 2 



f=u,d J A 



pUk F ,m 3 )^, (33) 



(34) 



where 



Ps( k U m f) = 1 



m f 



k<k' 



k 2 



(35) 



are the quark scalar densities. By imposing chemical 
equilibrium, supplemented by the conditions of charge 
neutrality and baryon number conservation, Eqs. 1241) - 
(|26|l along with Eqs. I|33|) - (|35|l form a system of six cou- 
pled equations, which are solved self-consistently to get 
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FIG. 1: Energy per baryon for symmetric matter calculated 
for the purely nucleonic case (solid line), with the MIT bag 
models (dotted and dot-dashed lines), and with the CDM 
(dashed line). 



X, c, k F ,k F ,k F , and k F . Once they are solved, we obtain 
the population of each quark flavor p/ = 7(fc£) 3 /67r 2 , as 
well as the one of the electrons. 



IV. RESULTS AND DISCUSSION 
A. Phase transition in symmetric matter 

Before studying the quark phase inside neutron stars, 
let us first discuss the EOS of symmetric matter. We 
calculate the EOS for cold symmetric nuclear matter in 
the BHF formalism with two-body and three-body forces, 
as described above in section [H] Then we calculate the 
EOS for u and d quark matter, using respectively Eq. 122|) 
for the MIT bag model and Eq. O for the CDM. The 
results are displayed in Fig. ^ The solid line represents 
the purely nucleonic case within the Brueckner many- 
body approach, whereas the different broken lines denote 
the results for the various quark matter EOS. We find 
that, in the range of baryon densities explored, the quark 
matter energy is always higher than that of symmetric 
nuclear matter, independently of the model adopted for 
describing the quark phase. Therefore nuclear matter 
is the favorite state. These results have been obtained 
choosing a s = 0, but we have checked that they are not 
very sensitive to the value of a s . 

As far as the results obtained with the MIT bag model 
are concerned, some clarifications are needed. In fact, 
apart from a bag constant B — 90 MeVfm -3 , we have 
also used a density dependent bag parameter B(p), whose 
parametrization is the one of Eq. (|28l) . This was adopted 
in previous works |2^|, and we do not repeat here the 
motivations and the technical details. 

The conclusion that indeed no transition to quark mat- 
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ter is present below (3, ...,4)po can be considered well 
established. At higher baryonic density the precise trend 
of the EOS of symmetric hadronic matter is surely more 
uncertain. A refinement of the EOS could lead to a tran- 
sition to a deconfined phase at a higher baryonic density, 
compatible with the value extracted from the CERN SPS 
experimental data (p c w 6po 1/fm 3 ). In particular, dif- 
ferent sets of three-body forces and the inclusion of three 
hole-line diagrams can produce indeed a stiffening of the 
hadronic EOS and the transition to quark matter can 
occur at this high density. A critical discussion on this 
point will be published elsewhere. However, these refine- 
ments do not affect the results in NS matter, since there, 
due to the appearence of hyperons and quark matter, 
the purely hadronic phase never reaches a density higher 
than (3, . . . ,4)po, and the uncertainity in the high den- 
sity EOS of nuclear matter plays no role. 

In order to facilitate the comparison between the dif- 
ferent quark models adopted in this work, it is useful to 
introduce in the CDM an effective bag parameter, which 
describes the difference between the quark matter energy 
density and the energy density of a system of free quarks 
of fixed mass m't having the same density and composi- 
tion. We thus define 



B e «(p) = ecDM - 7 



E E 

f=u,d,s k<k f 



k 1 



(36) 



where we take m' u = m' d = 0, m' s = 150 MeV in or- 
der to compare with the bag parameter used in the MIT 
model. As an alternative |52(, one can also define a bag 
parameter as the "non-quark" contribution to the energy 
density, 



B 'cs(p) = ecDM-7 E E V fe2 + TO /( cr 'X) 2 ( 3? ) 

f=u,d,s k<k f F 

= V( X ) + U(a,ir = 0). (38) 

The effective bag parameters B e g and B' eS are plotted 
versus the baryon density p in Fig. for symmetric mat- 
ter. We observe that, although at low baryon density the 
bag parameters used with the MIT model are larger than 
the one calculated with the CDM, asymptotically they 
all reach values in the range (50, . . . , 120) MeV fm -3 . 
However, in the former case the effective bag constant 
is a monotically decreasing function of the density, at 
variance with the CDM and also with the Nambu-Jona- 
Lasinio model, widely studied in Refs. jU EH • 



B. Phase transition in asymmetric /3-stable matter 

We now consider the hadron-quark phase transition in 
neutron stars, and calculate the EOS of a conventional 
neutron star as composed of a chemically equilibrated 
and charge neutral mixture of nucleons, hyperons, and 
leptons. The result is shown by the solid line in Fig.[3fa). 
The dashed line represents the EOS of beta-stable and 
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FIG. 2: The bag parameters B are shown as function of the 
baryon density for symmetric matter. The dot-dashed line 
shows the parametrization Eq. 1281 used in the MIT model, 
whereas the dashed lines represent the effective values of the 
CDM, calculated using Eqs. 136H or Eq. 138L respectively. 



charge neutral (u,d,s) quark matter obtained within the 
CDM, and the dotted and dot-dashed lines the results 
due to the MIT model with constant and density de- 
pendent B, respectively. Markers indicate the crossing 
points between the hadron and the quark phases. They 
lie inside the mixed phase region, whose range, at this 
stage, cannot be derived from the behavior of the energy 
alone. In spite of that, some qualitative considerations 
can be done. In particular, we notice that the phase tran- 
sition from hadronic to quark matter occurs at very low 
baryonic density when the CDM is used to describe the 
quark phase, whereas higher values of the transition den- 
sity are predicted with the MIT bag model. In fact, the 
density dependent bag parameter was introduced in or- 
der to shift this transition to high density and to explore 
the implications for the NS observables. 

A more realistic model for the phase transition between 
baryonic and quark phase inside the star is the Glenden- 
ning construction [2J, l54| , which determines the range of 
baryon density where both phases coexist, yielding an 
EOS containing a pure hadron phase, a mixed phase, 
and a pure quark matter region. In our previous papers 
|23| , we have used the Glendenning construction, demon- 
strating that in particular the influence on the maximum 
mass value is rather small. Apart from that, the realiza- 
tion of the mixed phase depends on the nuclear surface 
tension, which is currently an unknown parameter. 

Therefore, in the present work we adopt the simpler 
Maxwell construction. For that, we construct the phase 
transition from Fig.^b), showing the pressure as a func- 
tion of the baryonic chemical potential [1b- The coexis- 
tence region is determined by the intersection points be- 
tween the hadronic (solid line) and the different quark 
phases. We notice that the phase transition occurs at 
very low values of the baryon density, p w 0.05 fm -3 , 
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FIG. 3: Panel (a) displays the energy per particle vs. baryon density for beta-stable matter in the BHF approach (solid line) 
and for u,d,s quark matter obtained within the MIT model (dotted and dot-dashed lines) and the CDM (dashed line). Panel 
(b) shows the pressure as function of the baryonic chemical potential for all cases. 



when the CDM is used. The transition density that 
we obtain is somewhat lower than the value obtained in 
Ref. 0i where, using the Walecka model to describe the 
hadronic EOS, the mixed phase was found to occur in the 
interval 0.1 fm -3 < p < 0.31 fm -3 . However, the general 
result that the pure quark phase starts at a rather low 
density is in agreement with Ref. |47j . 

The phase transition constructed with the CDM is 
quite different from the ones obtained using the MIT 
bag model. In the latter case, the coexistence region 
is shifted to higher baryonic density, as seen in Fig.|3Jb). 
This implies a large difference in the structure of neutron 
stars. In fact, whereas stars built with the CDM have at 
most a mixed phase at low density and a pure quark core 
at higher density, the ones obtained with the MIT bag 
model contain a hadron phase, followed by a mixed phase 
and a pure quark interior. The scenario i s ag ain differ- 
ent within the Nambu- Jona-Lasinio model |53j , where at 
most a mixed phase is present, but no pure quark phase. 



Neutron star structure 



We assume that a neutron star is a spherically sym- 
metric distribution of mass in hydrostatic equilibrium. 
The equilibrium configurations are obtained by solving 
the Tolman-Oppenheimer-Volkoff (TOV) equations Q 



for the pressure P and the enclosed mass m, 
dP(r) Gm(r)e(r) 



dr 



[1 + P(r) /e(r)] [l + 47rr 3 P(r) /m(r)] 



1 - 2Grn(r)/r 



, (39) 



47rr 2 e(r) , 



(40) 



dm(r) 
dr 

being G the gravitational constant. Starting with a cen- 
tral mass density e(r = 0) = e c , we integrate out until 
the pressure on the surface equals the one corresponding 
to the density of iron. This gives the stellar radius R and 
the gravitational mass is then 



M G = m(R) = 4tt 



dr r 2 e( 



r) 



(41) 



We have used as input the equations of state discussed 
above for the CDM and the MIT bag model for the beta- 
stable quark phase, and the BHF for the hadronic matter. 
For the description of the NS crust, we have joined the 
hadronic EOS with the ones by Negele and Vautherin [5{| 
in the medium-density regime, and the ones by Feynman- 
Metropolis- Teller j5|| and Baym-Pethick-Sutherland [57| 
for the outer crust. 

The results are plotted in Fig. 0] where we display the 
gravitational mass Mq (in units of the solar mass M@) as 
a function of the radius R (left panel) and central baryon 
density p c (right panel). The dashed lines represent the 
calculation for beta-stable quark matter with the CDM, 
whereas the dotted and dot-dashed lines denote the re- 
sults obtained with the MIT bag model. Due to the use 
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FIG. 4: The mass (in units of solar mass Mq = 1.98 x 10 33 g) is displayed as function of radius (left panel) and central density 
(right panel), using different EOS. See text for details. 



of the Maxwell construction, the curves are not continu- 
ous 2]: For very small central densities (large radii, small 
masses) the stars are purely hadronic. Then a sudden in- 
crease of the central density is required in order to start 
the quark phase in the center of the star, corresponding 
to the phase diagram Fig.|2fb). Due to this effect, also a 
discontinuity in the possible radii arises. By performing 
the Glendenning construction, the curves would become 
continuous. 

We observe that the values of the maximum mass de- 
pend only slighty on the EOS chosen for describing quark 
matter, and lie between 1.5 and 1.6 solar masses. A clear 
difference between the two models exists as far as the 
radius is concerned. Hybrid stars built with the CDM 
are characterized by a larger radius and a smaller central 
density, whereas hybrid stars constructed with the MIT 
bag model are more compact, since they contain quark 
matter of higher density. 

This is also illustrated in Fig. |5] showing the differ- 
ent internal structure of stars with the MIT model and 
the CDM by comparing the populations of quarks and 
baryons inside a M = 1.4M© neutron star. At this value 
of the mass, within the MIT model only a thin hadronic 
layer is present, mainly composed of neutrons, followed 
by a small portion of crust, whereas in the case of the 
CDM one finds only crustal matter. The central (quark) 
density is substantially larger with the MIT model, in 
agreement with Fig.0J In both models, the abrupt tran- 
sition from quark to hadronic matter is a consequence of 
the Maxwell construction. 



V. CONCLUSIONS 

In this article we determine the structure of neutron 
stars, combining the most recent microscopic baryonic 
EOS in the BHF approach involving three-body forces 
and hyperons with different effective models describing 
the quark matter phase. 

Without allowing for the presence of quark matter, 
the maximum neutron star mass remains below 1.3 solar 
masses, due to the strong softening effect of the hyper- 
ons on the EOS, compensating the repulsive character of 
nucleonic TBF at high density. The presence of quark 
matter inside the star is required in order to reach larger 
maximum masses. 

We introduced a density dependent bag parameter 
B(p) in the MIT model in order to explore the maxi- 
mum NS mass that can be reached in this approach. We 
compare with calculations using a fixed bag constant and 
using the color dielectric model. Joining the correspond- 
ing EOS with the baryonic one, all three quark models 
yield maximum masses in the range (1.5,. ..,1.6) M@, 
while predicting slightly different radii. 

Our results for the maximum masses are in line with 
other recent calculations of neutron star properties em- 
ploying various phenomenological relativistic mean field 
nuclear EOS together with either effective mass bag 
model jH3 or Nambu-Jona-Lasinio model [5^ EOS for 
quark matter. 

The value of the maximum mass of neutron stars ob- 
tained according to our analysis appears rather robust 
with respect to the uncertainties of the nuclear and the 
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quark matter EOS. Therefore, the experimental obser- 
vation of a very heavy (M > 1. 7 M^ ) neutron star, as 
claimed recently by some groups [SJj (M ~ 2.2 Mq), if 
confirmed, would suggest that either serious problems are 
present for the current theoretical modelling of the high- 
density phase of nuclear matter, or that the assumptions 
about the phase transition between hadron and quark 
phase are substantially wrong. In both cases, one can 
expect a well defined hint on the high density nuclear 
matter EOS. 



FIG. 5: The particle populations inside a M = 1.4 Mq neu- 
tron star with the CDM (upper panel) and the bag model 
(lower panel) quark matter EOS. The markers indicate the 
beginning of the crust. 
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